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• • . Abstract 
> ■ 

^ I In the QCD Sum Rule determination of using the two-point correlator 

' of divergences of AS = 1 vector currents, the final uncertainty on nig is 

mainly due to the hadronic spectral function. Using a specific parameteri- 
zation which fully takes into account the available experimental data on the 
Ktc {1 = 1/2, = 0^) system, characterized by the presence of a relevant 
nonresonant component in addition to the resonant one, we find rns{l GeV) > 
120 MeV. In particular, varying only the parameters describing the nonreso- 
nant Kti component and we obtain fns{l GeV) = 125 — 160 MeV. This 
result is smaller than analogous ones obtained by using a parameterization in 
terms of only resonant states. We discuss how to systematically improve the 
determination of nig by this method. 



Light 'current' quark masses have an important role in the theoretical description of 
low energy hadronic physics, and their actual values are needed as an input to quantita- 
tively predict several relevant effects, such as, e.g., the breaking of chiral and vector flavour 
symmetries in hadron masses and transition amplitudes, and the size of CP violation in 
Kaon decays. 

Chiral perturbation theory represents a consistent framework where the (scale-inde- 
pendent) ratios of 'current' quark masses can be evaluated from experimental data. At 
the next to leading order in the chiral expansion one obtains [||, 0|: 

TTl Tfl 

— = 18.9 ±0.8 ; — = 0.553 ±0.043 . (1) 

However, for the determination of the individual light quark masses from the experi- 
mental data, one has to adopt some alternative non-perturbative method, such as lattice 
QCD or QCD sum rules. Once (at least) one of these masses is determined from these 
methods, the ratios in Eq. (|lD give access to the other ones. Here, we shall consider 
the determination of the strange quark mass from QCD sum rules. The reason for 
concentrating on (rather than on the other masses) is that there is more experimental 
information concerning the hadronic spectral function on which the QCD sum rule method 
is based, so that in this case one could expect the final result to be less model-dependent 
and more reliable. 

The present situation concerning the evaluation of is the following. Lattice QCD 
calculations, in quenched approximation, give for the running mass: 

msi2GeV) = 128±18Me\/ [3] (2) 
m,(2 GeV) = 100 ± 21 ± 10 MeV [| . (3) 

For = 380 MeV these values correspond to m«(l GeV) = 172 ± 24 MeV, and 

m,(l GeV) = 140 ± 29 ± 14 MeV, respectively. 

Regarding QCD sum rules, substantial progress has been achieved quite recently by 
the calculation to 0{a^) of the perturbative part of the two-point correlator of AS = 1 
scalar quark currents relevant to the determination of 0. Accounting for such 0{al) 
corrections, the value previously obtained in by using the 0{al) expression of the 
perturbative contribution to the sum rule:[| 

m,(lGe\/) = 189±32Me\/ [6] (4) 
^ Earlier QCD sum rules estimates, using different versions of the method, can be found, e.g., in refs. 
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has changed to the new value: 



m,(l GeV) = 203.5 ± 20 MeV [5] . (5) 

The uncertainties in Eqs. are related to the dependence of the result on the value 

of Ajj^^, which was chosen in the range 280 — 480 MeV , and to the variation of other 
input parameters in the sum rule analysis. In particular, as concluded in ||^, the largely 
dominant uncertainty is that on the hadronic spectral function, which in both ||^ and |^ 
is assumed to behave as the sum of two / = 1/2, J-'^ = 0^ (Kit) resonant states. 

Following this observation, we would like to discuss a parameterization of the relevant 
hadronic spectral function which, at least in principle, should be able to fully exploit 
the current experimental information on the scalar J = 1/2 channel. In particular, this 
parameterization takes into account the fact that, in addition to the resonance component, 
there exists also a nonresonant continuum. As typical of s-wave channels, the nonresonant 
component persists rather far from the threshold energy range and can interfere with the 
resonances Our example shows that the inclusion of the nonresonant contribution 
may give a significant effect on the sum rule and, also, indicates some directions in order 
to improve the determination of in this approach. 

For convenience, we briefly describe the main points of the QCD sum rule derivation 
of rris we are interested in. The basic quantity is the two-point correlator 

^(g2) = ijdx e^«^ < 0\T[J{x)j\0)]\0 > (6) 

where J{x) = i{ms —mu)s{x)u{x) is the divergence of the AS* = 1 vector current V^{x) = 
s{x)'y^u{x). The second derivative {q'^) = {d"^ / {dq^Y)'^ [q^] obeys an unsubtracted 
dispersion relation: 

*V)=2rd,^4fi— , (7) 

Jo \s — q — "^ej 
with the spectral function p(s) given by 

p{s) = -Im ^(s) . (8) 

TT 

Particularly useful for low-energy phenomenology is the Borel transform ^f" [M"^) ||11|| , 

defined by the application of the operator B{M^) = ; — (— — --)"" in the limit n ^ oo, 

n! ^d[—q^) 

—q^ — >■ oo, —q^/n = = const to Eq. (^. This results into: 

V{M^) = ^l^dsp{s)e-^/'^' . (9) 
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The main advantage of this transformation is that, due to the presence of the exponential, 
for moderate values of the Borel parameter (of the order of one to a few GeV'^) the 
r.h.s of Eq. @ is mostly sensitive to the spectral function p{s) in the low energy range, 
where it can be computed in terms of the available experimental information. 

The quark 'current' mass enters into the l.h.s of Eq. (|^) as a parameter, when the 
correlation function is evaluated in QCD for = —q^ ^ ^qcd t>y means of 

the Operator Product Expansion. Once p{s) is parameterized by a hadronic representa- 
tion inferred from experimental data, the resulting Eq. (P) represents a (M^-dependent) 
equation relating iris to experimental data and known QCD parameters. 

To arrive at the final result, one applies the notion of global quark-hadron duality, 
which essentially consists in identifying dsp{s)e~'^/'^^^ = dspopE{s)e~'^^^^^ , since, 
due to asymptotic freedom, quark and gluon degrees of freedom (rather than hadrons) 
should dominate above an effective energy threshold Sq. Accordingly, Eq. (P) takes the 
form 

^"iM')\oPE = j^l ds p{s) e-^"^' + — j^^^ dspoPEis) e-'l^'^ , (10) 

and the ultimate numerical determination of will be the one for which this relation is 
stable in the, as yet undetermined, parameters and Sq. 

The OPE expression for \I^"((5^) can be given in terms of a perturbative and a non 
perturbative contribution: 

*"(Q') = ^;(Q')+*;p(Q') • (11) 

For illustrative purposes, we report here just the leading order expression of ^"pi^Q"^): 
3 (m^(/i) - m„(yu))2 f \\as{p) „«s(/i) , 



^ Stt^ 3 tt tt p? ' 

6 ms^{p){^s{[i) -mu{p)f ^ ^ 28 a,(/i) ^ ^^M ipg*^^] ^2) 

and we refer to |^ for the explicit, lengthy expressions of the order and a\ contribu- 
tions.0 In Eq. (|T^), /i is an a priori arbitrary renormalization mass scale. Since \l/"((5^) is 
related to a physical observable, it obeys the homogeneous renormalization group equa- 
tion: 

/^|^*"(Q') = , (13) 



We keep niu 7^ only in the coefficient of the perturbative function 
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and therefore the scale dependence of the renormahzed parameters as and appearing 
in the perturbative calculation of \E'"(Q^) must cancel against log/i factors also appearing 
in \E'"(Q^). The Borel transform of \E'"((5^) in the approximation of Eq. (|12D is given by: 

^„ _ 3 (m.(ri-m„(^))^ ^ ^ llaM _ j^Uff) (log^^ + ^(l))} 
_^ m,\^)(mM -mJp)f .^ 28^ _ 4£iM(iog + .^.(2))} .(14) 

where is the dilogarithmic function. The choice fi = M allows to resum the loga- 
rithmic terms, transforming the dependence on /i of the running mass m(/i) and of the 
running coupling constant ^^(/i) into a dependence on the Borel parameter M. 

In an analogous way one can write the non perturbative contribution to 
Referring to for the detailed expression of the operator expansion, we report here only 
the Borel transformed contribution of D = 4 operators: 

^7Vp(M') = M6 {2<m,n^>o (l + ^(y-2^(l)-21og— )) 

3 _4/7r 5 15 15, 

/s and Ig are RG invariant combination given by (for rtf = 3): 

/, = m, < ss >o +-^m^(— - ^) (16) 

+ 7^(1 + • (17) 

The hadronic contribution to the spectral function p{s) can be obtained by inserting 
a set of intermediate states with = 0^ and I = j into the correlator (^). The 
simplest examples are the two-particle states \K7c >, \Kr] >, \Kri' >. In particular, the 
contribution of the {Kn > intermediate state, which is expected to be the dominant one 
and whose features are better known from the theoretical as well as the experimental 
point of view, f\ can be written as: 



(Kn), ^ 3 \/(s-s+)(s-s_) 

OZTT^ S 



^Multiparticle states should be suppressed by phase space. 
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where s± = {Mk ± Mn-)^, and d{s) is the scalar form factor related to the K^s decay form 
factors f±: 



1 



< n\p')\sj,u\K+{p) >= [(p + p')J^ + {p-p')J_] 



V2 



d{s) = {Ml-Ml) U{s) + 



Ml-Ml- 



Us) ={Ml-Ml)Us) . 



(19) 



(20) 



From the theoretical point of view, d{s) and fo{s) can be considered as analytic functions 
of the complex variable s, with a cut on the real axis starting at the threshold s+ = 

Furthermore, in the range < s < s_, which is the physical one for K^^ decay, /o(s) 
admits a linear expansion for small s: 

/o(.) = /o(0)(l + Ao^) (21) 

where, from one-loop chiral perturbation theory 0, |T2|, one has the predictions: 

/o(0) = 0.973; Aq = 0.017 ± 0.004 . (22) 



The theoretical value of the slope Aq in Eq. (|22D agrees with the experimental result from 
the high statistics analysis of K^^ decays, which gives the result Aq = 0.019 ± 0.004 ||T3[|.P| 
The linear extrapolation of Eq. ( pi]) with the theoretical values (|2^) from the decay region 
to the threshold s+ would imply the prediction d{s^) = 0.30 ± 0.02. 



As a final constraint, inspired by the quark counting rules ||I5[, we assume the asymp- 
totic s-behaviour: d{s) ~ 1/s. 

Experimental information on the Ktt system was obtained from the analysis of the 
reaction K~p — > K^ir^n some time ago [|16|. The partial wave analysis of Kit — > Kit 
provides evidence, in the 0"*", 1 = 1/2 channel of a well-established i^Q(1430) resonance, 
with Mr = 1429 ± 4 ± 5 MeV and Tr = 287 ± 10 ± 21 MeV, and a signal for a not yet 
confirmed i\:*(1950) state, with Mr = 1945 ± 10 ± 20 MeV and = 201 ± 34 ± 79 MeV 



161 , [T^. Moreover, a non negligible nonresonant component underlying the ^^0(1430) 



resonance shows up in the measured low energy Kit phase shifts. 



^ However, this value seems not quite in agreement with the resuh from K^^ decays p^ . 
^ Actually, a small curvature from higher order terms in the low energy representation of /o is admitted, 
and in principle might be perceptible at s+ which is much larger than the decay endpoint s_. 
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In 0, 0, the behaviour of |(i(s)p appearing in Eq. (|TB|) is modeled by the sum of two 
Breit-Wigner forms (with masses and widths as reported previously), normalized at s = s+ 
to the theoretical prediction from (|2l|) and (0). To model an alternative parameterization 
of (|T^) which includes all the experimental information mentioned above, in particular 
the existence of the nonresonant component, one can attempt a construction of the form 
factor d{s) based on analyticity properties and asymptotic behaviour, with the measured 
Kir phase shifts as an input, consistently with the final state interaction theorem. Such 



a construction of d{s) can be realized by assuming the following representation ||T8 



with d{0) determined from Eqs. (pO]), (p^) and 6{s) the Kn (/ = ^, = 0+) phase shift. 
As usual in the applications of this representation, we do not consider the possibility of 
zeroes for d{s), which would require also a polynomial factor in (0). 

The Kit phase shift is well known in the range of invariant mass from s+, to (1.7 GeV)'^: 
it can be parameterized as the sum of an effective range formula and a resonant phase: 

S{s) = 6er{s) + 6bw{s) (24) 

with: 

Ser{s) = arctg [aq{s){l + bq\s))] (25) 
(g(s) is the Kit CM momentum) and 

MnTnis) 

Obw = (^rctg[—^ J , (26) 

— S 

where Mji is the mass and Tji{s) is the s-dependent width of the A'q(1430) state: 



Mr q{s) 
rsq{Ml)- 



^nis) = ^^Tn (27) 



The parameters a and h have been fitted, with the result [Q: 

a = 2MGeV-^ , b = -1.37 GeV~'^ . (28) 

As for the region s > (1.7 GeVy, above the Krj' threshold, inelastic effects are observed 
161 , |10| .0 The inclusion of such effects would require a coupled channel analysis of the Kn 



^The contribution of the Krj state is flavour-S'C/ (3) suppressed ||lC 
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and Kt]' states, with further contributions to p^^^^\ Since we do not consider here such 
additional contributions, strictly speaking the result of our analysis is a lower bound for 
rUs, due to the positivity properties of the spectral function p(s). However, the exponential 
factor in Eq. (j^) should suppress the contribution of higher states. This is confirmed by 
the numerical analysis of ref. 0, where the contribution of the i^'Q(1950) has a very small 
influence on the result for m^. 

As far as d{s) is concerned, the asymptotic 1/s behaviour at large s can be obtained 
from if 6{s) 180°, and in this regard we fix 6{s) = 180° for s > (1.7 GeV)^. This 
is a delicate assumption from the numerical point of view, which, however, is supported 



by the available experimental data |T0|, [16 



We can remark that, using the parameterization (p^D of the phase shift, the form 
factor d{s) in (p3|) reproduces the slope predicted by chiral perturbation theory, as already 
observed in and shown in Fig. 1. This feature is not obtained, in the framework of the 
representation ( PB] ) for d{s), if the phase shift is parameterized in terms of the resonant 
^Bw{s) phase only. Therefore, in the approach considered here, the inclusion of the 
nonresonant component is needed on phenomenological grounds. The obtained spectral 
function is depicted in Fig. 2, where we compare the result from the parameterization 
(IH)"(H) with the case of the pure Breit-Wigner form obtained using a = b = 0. f\ The 
substantial reduction of the resonance peak shows that even a moderate nonresonant 
contribution in the Kn phase shift 6{s) can generate a significant variation of the spectral 
function via the exponential form of the analytic representation (pSj). 

At this point, we perform the numerical analysis of the sum rule, following the same 
procedure adopted in |^. The only slight difference with respect to |^ is that, in the 
expansion of the f3 function relevant to the perturbative part of the sum rule: 



/3(a) = ^/?„a" (« = -) (29) 



we use the coefficient (34 recently computed in the MS scheme [^, namely the set of 
values: 

9 ^ „ ^ 3863 ^ 281198 890^,^, 

Pi = — , P2 = -8, p3 = , p4 = C 3) (30) 

PI 2' '^^ ' 192 ' 4608 32 ^ ^ ' 

for Nc = 3 and Uf = 3; ( is the Riemann zeta function. The numerical value for jS^ differs 

by a factor of two from that used in ||^ obtained by a Fade approximant |21]. 



^The identity of the representation (|2^) with a Breit-Wigner form in the case of a pure resonance can 
be shown directly [p^ . 
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Moreover, for the analogous expansion of the anomalous dimension 7(a): 7(a) = 
SnTnfl", we use the first three computed coefficients (for Nc = 3, Uf = 3) 

91 8885 

7i = 2, 72 = — , 73 = - 5C(3) , (31) 

whereas for 74 we use the result of two Pade approximants 

7(a) ^ A + Ba 
a 1 + Qa ^ ^ 

{A, B, Q numerical coefficients) giving 74 = 73/72 f^, |22|, and 

a 1 + Ka + ba^^ 

giving 74 = 27273/71 — 72/71- The difference between the two approximants (|3^) and 



(|33|) is about 3%, and has no practical influence on the final result for m^. As for the 

parameters in the theoretical side of the sum rule, we use the values: < uu > |^=i ceV = 

(-0.225 GeVf, < -G^ >= (2 - 6) 10"^ GeV^ and — = 0.7 - 1.0; in the coefficient 

vr <uu> 
of eq. (|12D we use = ms/34.2 as from Eq. ([l|). The dependence of the result for 

from these nonperturbative contributions is rather weak, and therefore the uncertainty 

corresponding to the variation of the parameters in the considered ranges is found to be 

small, of the order of 1 MeV. 

We varied the threshold Sq in the range 5 — 7 GeV"^, the Borel parameter in the 
range 1 — 9 GeV"^, and we considered the three values for A^(^^: 280, 380 and 480 MeV. 
Specifically, the best stability in of the resulting is obtained for sq = 5.5 GeV^ 
when A^^ = 380 MeV; for sq = 5.0 GeV^ when A^^ = 280 MeV and for sq = 6 GeV^ 
when A^^ = 480 MeV. f\ 

The stability curves for the running mass ms(l GeV) are depicted in Fig. 3 where only 
the dependence on A^^^ and is displayed. Changing only A^f^^ the result for 
changes in the range 130 — 140 MeV. An additional uncertainty is due to the parameters 
describing the continuum Kn component; changing the parameter a in the effective range 
formula (|25|) by 20%, and the other parameters in the considered ranges, we have 

m,(l GeV) = 125 - 160 MeV (34) 



^After this paper was completed, a calculation of 74 in the MS scheme appeared in . The computed 
value: 74 = 88.5258 differs from the results in ( |3^ , ^3|) by 12%. This effect has no influence on the result 
for nis. 

^ Concerning the possible instanton contribution to the correlator (^) , even though a quantitative 
assessment is difficult, we argue that such contribution is within the uncertainty related to the hadronic 
spectral function, mainly due to the large energy scale involved in the calculation of the correlator Eq. (^. 
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with a lower bound 

m,(l GeV) > 120 MeV. (35) 

The value (0) derives from a parameterization of the hadronic spectral function which 
has the correct analyticity properties and uses all the available experimental information. 
This determination can be systematically improved by a dedicated analysis of the scalar 
/ = I Kn channel using, e.g., the semileptonic decays of the r lepton. 

In principle, the result ( |34D is reflected in the lower bound (|35|) due to the neglect 
of higher states; however, we do not expect that the actual determination of in the 
theoretical framework considered here can significantly differ from it. The value (0), 
lower that the result obtained assuming the resonance dominance, shows the significancy 
of the Kit nonresonant continuum. Q 



^"Similar conclusions have also been reached in ref. |25| 
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FIGURE CAPTIONS 



Fig. 1 

Form factor d{s) obtained using the representation Eq. (|23|) (continuous line). The xPT 
prediction Eq. (^1]), with Aq = 1.69 10~^, is also shown (dashed hne). 



Fig. 2 

The hadronic spectral function p(s) obtained including only the resonant i^Q(1430) state 
(upper dashed line), or using the representation Eq. ( p3D (lower lines); the three low-lying 
curves correspond to a 20% variation of the parameter a in Eq. (p5|) . 



Fig. 3 

The running mass rna{^) at the scale /i = 1 GeV as a function of the Borel parameter 
M2, for = 380 MeV (continuous hne), = 280 MeV (dotted line) and 

iCff = 480 MeV (dashed line). 
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